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1.  INTRODUCTION 


The  Air  Force  Strategic  High-Altitude  Radiance  Code  (SHARC)^  has  been 
developed  to  address  the  strategic  requirements  for  modeling  infrared  (IR)  background 
radiation  in  the  upper  atmosphere.  One  of  the  most  critical  applications  is  in  modeling 
spatial  variability,  which  is  an  important  effect  in  the  solar  terminator  (earth  shadow  edge) 
region.  Two  types  of  effects  can  occur  which  lead  to  spatial  radiance  variations  at  the 
shadow  edge.  One  is  a  change  in  molecular  number  densities  due  to  photochemical 
processes,  such  as  the  photolysis  of  ozone  (O3).  The  other  is  a  change  in  the  vibrational 
excitation  with  the  solar  zenith  angle.  The  net  result  in  either  case  can  be  a  substantial 
variation  in  the  spectral  radiance  for  many  molecular  bands,  including  O3  (9.6  ^m),  CO2 
(4.3  fim  and  2.7  /tm),  and  OH  (~  1.9  ^m  and  ~3.3  /xm). 

The  most  recent  version  of  SHARC,  SHARC-3,  permits  accurate  modeling  of  the 
terminator  region  through  the  use  of  multiple  atmospheric  profiles  corresf>o,iding  to 
different  solar  angles.  The  spectral  radiance  may  be  calculated  for  an  arbitrary  line  of 
sight  (LOS)  passing  through  up  to  seven  profiles.  Interpolation  is  used  to  characterize  the 
atmosphere  between  the  solar  angles  assigned  to  each  profile.  The  number  of  profiles  is 
sufficient  to  accurately  describe  the  characteristics  of  the  shadow  edge  in  any  arbitrary 
view  (parallel,  perpendicular,  or  oblique). 

In  the  current  investigation,  SHARC-3  was  used  to  simulate  terminator  limb 
radiances  for  CO2,  O3,  and  OH  that  correspond  to  field  observations.  Dawn  terminator 
data  for  the  CO2  4.3  /xm  band  were  selected  from  the  1977  SPIRE  rocket  measurements.^ 
New  SHARC  simulations  were  performed  using  an  input  atmosphere  from  the  SHARC 
Atmospheric  Generator  (SAG),  a  model  that  utilizes  recent  climatology  databases, 
empirical  models,  and  theoretical  calculations  to  construct  appropriate  atmospheric  profiles 
for  any  desired  conditions.  The  calculated  CO2  radiances  are  in  good  agreement  with  the 
SPIRE  data  as  well  as  with  simulations  performed  using  GL’s  line-by-line  code  RAD.^  In 
addition,  we  have  simulated  dusk  terminator  emissions  from  O3  and  OH  measured  in  the 
1991  CIRRIS-IA  space  shuttle  experiment.^  The  reasonable  agreement  obtained  serves  to 
validate  both  the  SHARC  code  and  the  SAG  solar  angle-dependent  model  for  O3 
concentration. 


2.  CALCULATION  METHODS 


Two  slightly  different  methods  for  performing  SHARC  terminator  calculations  were 
used  in  the  current  investigation.  The  CO2  4.3  /xm  simulations  for  SPIRE  were  performed 
using  a  multiple-regions  version  of  SHARC-3.  Five  atmospheric  regions  corresponding  to 
different  solar  angles  were  defined  with  explicit  geographic  boundaries.  For  convenience, 
a  transformed  coordinate  system  was  used  in  which  the  scan-center  tangent  point  was 
placed  on  the  equator  and  the  solar  angle  served  as  the  longitude  coordinate.  (See  Table 
1.)  The  atmospheric  regions  were  stacked  along  the  equator  and  bounded  by  longitude 
Fnes;  the  boundary  longitudes  were  chosen  as  the  mean  of  the  solar  angles  assigned  to  the 
two  regions  on  either  side  of  the  boundary.  One  of  the  five  regions  was  "extended",  i.e., 
had  no  boundary.  This  region  was  placed  at  one  end  of  the  stack. 

For  simplicity,  the  vehicle  was  regarded  as  stationary  during  each  spatial  scan  (this 
is  a  good  approximation  for  SPIRE  and  other  rocket  experiments).  Since  the  azimuth 
angle  wa>  also  fixed  during  the  scan,  the  entire  limb  radiance  profile  for  the  scan  could  be 
mapped  out  by  changing  only  the  LOS  tangent  height  in  the  SHARC  input  file.  This  was 
done  in  an  automated  fashion  by  running  SHARC  from  an  external  routine.  The  rout  ne 
modified  the  tangent  height  in  the  input  file,  executed  SHARC  to  compute  the  LOS 


Table  1.  SPIRE  Viewing  Coordinates  in  the  Transformed  System 
Used  for  the  Multiple-Regions  SHARC  Calculations. 


Vehicle 

Scan 

Height 

Longitude 

Latitude 

1 

188.57 

92.40 

-10.367 

4 

250.09 

92.08 

-9.922 

8 

279.00 

91.80 

-0.0675 

9 

266. 1 1 

91.24 

-2.5958 

10 

235.97 

91.03 

-4.9565 

11 

213.26 

90.93 

-7.2198 

LOS 

Azimuth 

Height 

Tangent  Point 

Longitude 

I^ititude 

26.08 

56.60 

97.43 

0 

39.16 

87.76 

100.07 

0 

89.65 

158.43 

102.73 

0 

-75.93 

153.63 

81.00 

0 

-62.13 

125.10 

81.75 

0 

-37.37 

131.18 

85.45 

0 
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spectral  radiance,  integrated  the  result  within  the  selected  bandpass,  and  wrote  an  output 
file  containing  the  bandpass  radiance  and  tangent  height  values. 

The  O3  and  OH  simulations  for  CIRRIS-IA  were  performed  using  the 
multiple-profiles  feature  in  the  July,  1993  release  of  SHARC-3.  Seven  atmosphere 
profiles,  corresponding  to  a  92°-104°  solar  zenith  angle  (SZA)  range,  were  specified  in  a 
single  quiescent  region.  The  actual  latitudes  and  LOS  azimuth  angles  for  the  experiment 
were  used.  The  correct  placement  of  the  sun  was  achieved  through  the  choice  of  input 
longitude. 

The  only  real  difference  between  using  a  single  profile  within  multiple  regions  (in 
the  CO2  calculations  for  SPIRE)  as  opposed  to  multiple  profiles  within  a  single  region  (in 
the  calculations  for  CIRRIS-lA)  is  that  in  the  latter  case  the  vibrational  state  populations 
are  interpolated  along  the  LOS,  whereas  in  the  former  case  they  are  constant  within  each 
region  and  change  abruptly  at  the  boundaries.  For  radiators  such  as  CO2  where  the  solar 
angle-dependence  of  vibrational  state  populations  is  gradual,  the  calculated  limb  radiances 
are  similar  with  either  approach.  With  O3  and  OH  emissions,  which  change  dramatically 
over  a  several-degree  solar  angle  range,  the  interpolation  provided  in  the  multiple-profiles 
method  should  give  smoother  results. 
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3.  RESULTS  AND  COMPARISONS 

3. 1  CO2  4.3  /xm 

Calculations  of  the  4.19-4.49  /xm  bandpass  radiance  were  carried  out  to  simulate  the 
SPIRE  "terminator/night"  scans  1,  4  and  8  and  the  "day"  scans  9  and  11.  For  the  former, 
five  layers  at  2°  SZA  increments  were  used,  for  a  total  SZA  range  of  10°  along  the  LOS. 
This  is  sufficient  for  coverage  of  tangent  heights  up  to  150  km  or  greater.  For  scan  8, 
v/hich  looked  much  farther  into  the  night  than  the  other  two,  a  low-solar-angle  extended 
region  was  substituted  for  the  daylit  extended  region  used  in  the  scan  1  and  scan  4 
calculations.  For  the  day  scans  9  and  11,  which  have  generally  greater  East-West  ranges, 
the  atmospheric  regions  were  made  three  and  four  degrees  wide. 

3.1.1  CO2  Vibrational  Temperatures 

To  obtain  some  insight  into  the  SZA-dependence  of  the  CO2  excitation  and  the 
convergence  properties  of  the  calculations,  we  tabulated  the  vibrational  temperature  of  the 
main  isotope  (00011)  state  (associated  with  the  cold  band)  for  each  atmospheric  region. 
As  expected,  there  is  a  general  decrease  with  increasing  SZA  (Figure  1).  The  temperatures 
for  one  set  of  altitudes,  around  119  km,  track  one  another  closely.  High  temperatures  in 
region  1  are  followed  by  approximately  level  temperatures  in  the  darker  regions. 

A  similar  trend  appears  at  lower  altitudes  around  91  km,  but  the  temperatures  do 
not  track.  The  small  fluctuations  appearing  at  this  altitude  are  statistical,  stemming  from 
the  finite  number  of  Monte  Carlo  "photons"  used  in  the  NEMESIS  routine  that  calculates 
the  vibrational  populations.  Ten  thousand  photons  per  vibrational  band  were  used  to 
develop  the  Figure  1  results,  and  should  be  adequate  for  most  applications.  If  this  number 
is  decreased  to  N  =  500,  much  larger  fluctuations  occur,  as  in  Figure  2.  The  relative 
magnitudes  of  the  fluctuations  in  these  two  figures  are  consistent  with  the  statistical  "v/n 
scaling,  as  expected. 
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Figure  2.  Same  as  in  Figure  1  But  Using  500  NEMESIS  Photons. 


3. 1 .2  Radiance  Profiles 


The  4.3  nm  radiance  profiles  from  SHARC  appear  in  Figure  3  along  with  the 
SPIRE  data.  Since  the  data  display  significant  scatter  not  correlated  with  SZA,  data  for 
similar  angles  have  been  collected  into  two  sets,  corresponding  to  day  and  night  sides  of 
the  terminator.  The  scan  1  data  are  very  limited  and  are  not  shown.  Two  SHARC  profiles 
are  shown  for  each  set  to  indicate  the  sensitivity  of  the  caU  ulations  to  the  SZA.  No 
differences  are  found  between  the  various  nighttime  calculations  below  100  km,  where 
strong  atmospheric  absorption  of  sunlight  leads  to  total  shadow.  In  other  cases,  minor 
differences  are  found.  The  agreement  between  the  SHARC  calculations  and  the  SPIRE 
data  is  seen  to  be  within  the  data  scatter.  As  can  be  seen  in  Fig.  3,  the  SHARC 
calculations  for  99  and  107  SZA  are  essentially  nighttime  except  at  the  higher  altitudes. 


O 


In-Band  Radiance  (W/cm^/sr) 


Figure  3.  Comparison  of  4.3  /xm  Limb  Radiance  Measurements  from  the  SPIRE 
Rocket  Experiment  with  SHARC-3  Simulations. 


Preliminary  SHARC  calculations  have  also  been  performed  to  simulate  CIRRIS  lA 
data,  which  have  better  signal-to-noise.  The  results,  which  will  be  reported  in  the  future, 
indicate  agreement  to  within  a  factor  of  2  over  the  50  to  125  km  tangent  height  range.  The 
SHARC  radiances  tend  to  be  lower,  particularly  at  the  higher  tangent  heights.  An 
underprediction  of  this  magnitude  has  also  been  found  by  Wintersteiner  and  co-workers  at 
nighttime  using  a  different  radiation  code  (RAD)^  with  similar  kinetic  mechanisms. 

3.2  O3  9.6  fim 

In  contrast  with  CO2,  the  variation  in  ozone  radiance  at  the  solar  terminator  is  due 
to  an  abrupt  change  in  concentration.  Thus,  the  SHARC  terminator  results  depend 
critically  on  how  the  ozone  concentration  profile  changes  with  solar  zenith  angle  in  the 
model  atmospheres.  A  discussion  of  the  ozone  profile  derivation,  exc.erpted  from  the 
SHARC  Atmosphere  Generator  (SAG)  documentation,^  is  presented  below. 

3.2.1  Derivation  of  Ozone  Concentration  Profile 

The  ozone  profile  is  based  on  a  combination  of  the  NRL  trace  species  climatology 
database^  and  a  photochemical  equilibrium  model,  described  below.  The  NRL  database 
contains  daytime  mixing  ratios  and  night/day  ratios,  and  is  used  below  72  km.  The 
latitude  range  is  restricted  to  ±60”  to  avoid  interpolation  difficulties  associated  with  the 
onset  of  continuous  polar  day  or  night.  The  photochemical  model  is  used  at  80  km  and 
above,  and  the  two  regions  are  connected  by  logarithmic  interpolation. 

The  photochemical  equilibrium  model  accounts  for  ozone  formation  from  O  atoms 
and  destruction  via  H  atom  and  O  atom  reactions  (the  former  yielding  OH)  and  solar 
photolysis.  The  inputs  are  the  photolysis  rate  and  the  temperature,  N2,  O2,  O,  and  H 
profiles,  which,  as  discussed  in  the  SAG  documentation,  are  based  almost  entirely  on  the 
MSISE-90  model. ^  Along  with  the  standard  three-body  recombination  reaction,  a 
two-body  radiative  association  term  is  also  included,  and  has  a  significant  and  increasing 
impact  above  100  km.  The  value  of  the  rate  constant,  8x10'^*  cm^/molec/s,  is  estimated 
from  rate  constants  for  formation  and  relaxation  of  the  O2-O  complex®  and  from  1^3  band 
Einstein  coefficients.  It  is  almost  certainly  an  upper  limit.  We  feel  that  it  is  safer  to  risk 
overestimating  the  ozone  than  to  underestimate  it.  Furthermore,  any  overestimation  from 
this  2-body  channel  would  be  reduced  by  the  presence  of  certain  62*+  O2  reactions  that 
have  been  proposed  to  form  ozone^  but  which  are  not  presently  included  in  SAG. 
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The  transition  between  day  and  night  profiles  is  handled  using  a  simple  empirical 
model  that  approximates  the  terminator  behavior  in  the  time-dependent  calculations  of 
Rodrigo  et  al.,*^  which  were  performed  for  mid-latitude  (45°)  equinox  conditions  above  60 
km.  The  transition  is  made  by  linearly  interpolating  over  a  5°  range  of  SZA  computed  for 
0.075  hr  earlier;  the  0.075  hr  (4.5  min)  difference  accounts  for  the  approximate 
photochemical  lag  time.  The  midpoint  of  the  interpolation  region  is  taken  to  vary  linearly 
with  altitude  at  the  rate  of  0. 12°  per  km. 

3.2.2  Ozone  Profile  Comparisons 

Figure  4  compares  the  daytime  and  nighttin  concentrations  from  SAG  and  from 
the  NRL  climatology.  The  SAG  and  NRL  climatology  results  differ  above  72  km  where 
the  transition  to  the  photochemical  model  is  made  in  SAG,  but  the  difference  is  generally  a 
factor  of  two  or  less.  According  to  the  photochemical  model,  above  80  km  ozone  is 
destroyed  almost  exclusively  by  reaction  with  H  atoms.  Much  of  the  difference  between 
the  NRL  and  SAG  nighttime  profiles  may  have  to  do  with  inaccuracy  in  the  MSISE-90  H 
atom  profile  used  by  SAG.  Recent  observations'*  suggest  that  the  MSISE-90  profile  may 
be  somewhat  low,  especially  near  80  km.  Further  investigation  of  this  issue  is  warranted. 

The  temporal  behavior  of  the  SAG  ozone  concentrations  in  the  dawn  and  dusk 
terminator  regions  is  shown  in  Figures  5a  and  5b.  Comparing  them  with  Rodrigo  et  al.’s 
calculations  in  Figures  6a  and  6b,  it  is  seen  that  while  the  absolute  concentrations  from 
SAG  are  different,  particularly  at  night,  both  the  day/night  transition  times  (indicating  the 
shadow  edge  location)  and  the  slopes  of  the  curves  (indicating  the  sharpness  of  the  edge) 
are  very  similar  throughout  the  60-100  km  altitude  range.  At  lower  altitudes,  both  the 
NRL  climatology  database  and  our  own  survey  of  literature  measurements  indicate  that, 
contrary  to  Rodrigo’s  calculations,  the  diurnal  variation  of  ozone  is  veiy  small.  Therefore, 
any  error  in  the  terminator  edge  model  will  be  negligible  below  60  km. 

In  conclusion,  the  SAG  model  should  provide  a  good  empirical  description  of  ozone 
concentration  in  the  daytime,  nighttime,  and  terminator  region. 
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Comparison  of  SAG  Model  and  NRL  Climatology  Ozone 
Concentrations  for  Fall  Equinox,  45°  N  Latitude,  Moderate-Activity 
Conditions. 
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LOCAL  TIME  (hrs) 


LOCAL  TIME  (hns) 


Figure  5.  Time- Dependent  Ozone  Concentrations  from  SAG  in  the  Solar 

Terminator  Regions  for  Fall  Equinox,  45*^  N  Latitude, 
Moderate-Activity  Conditions. 


Figure  6.  Time-Dependent  Ozone  Concentrations  from  Rodrigo  et  al.’°  in  the 
Solar  Terminator  Regions  for  Fall  Equinox,  45®  N  Latitude. 
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3.2.3  Terminator  O3  Radiance 

SHARC  calculations  were  performed  to  simulate  selected  data  from  the  CIRRIS-IA 
experiment  on  O3  radiance  in  the  dusk  terminator.’’  The  data  were  taken  from  a 
horizontal  scan  near  20®  S  latitude  during  Block  PCI  IG,  when  a  SZA  range  of  93®  to  101" 
was  covered  with  nearly  constant  elevation  and  azimuth  angles.  During  this  period,  two 
radiometer  detectors,  numbered  1-4  and  5-6,  acquired  data  at  89  km  and  81  km  tangent 
heights,  respectively,  in  the  8.0-11.9  /xi.i  bandpass  The  LOS  was  nearly  parallel  to  the 
terminator. 

The  data  and  SHARC  simulations  are  shown  in  Figure  7,  The  SHARC  simulations 
show  a  smaller  altitude  dependence  on  the  night  side  and  slightly  underestimate  the  day 
side  radiance,  but  on  the  whole  they  agree  reasonably  with  the  data.  The  agreement  is 
better  at  81  km,  which  is  closer  to  the  peak  of  the  ozone  emission,  than  at  89  km.  The 
differences  on  the  night  side  may  have  to  do  with  a  mismatch  between  the  calculated  and 
actual  ozone  peak  altitude,  which  averages  around  83  km  in  CIRRlS-lA  nighttime 
measurements  but  can  fluctuate  by  several  km.  For  these  particular  data,  the  peak  appears 
to  be  several  km  lower  than  average.  At  89  km,  the  .scale  height  of  the  emission  is  very 
steep,  so  that  a  small  change  in  the  peak  altitude  can  have  a  large  effect  on  the  radiance. 
Uncertainties  in  the  CIRRIS-IA  pointing  could  also  explain  some  or  most  of  the 
differences.  The  tangent  height  uncertainty  is  estimated  to  be  1  to  2  km  for  these  data. 
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3.3  0HAV=1 


The  OH  AV  =  1  Meinel  band  emission  above  80  km  is  generated  from  the  H  +  O3 
reaction.  The  H  and  O3  concentration  profiles  have  been  discussed  in  Section  3.2.1.  The 
OH  emission  is  expected  to  show  a  similar  terminator  behavior  as  the  O3  emission,  since 
the  H  atom  concentration  should  not  vary  across  the  terminator  at  these  altitudes.  Not 
surprisingly,  we  find  the  SHARC/CIRRIS-IA  data  comparison  for  OH  to  be  similar  to  that 
for  O3. 

3.3.1  Terminator  OH  Radiance 

Figure  8  shows  2. 5-3. 4  fim  bandpass  data  from  CIRRIS-IA  that  were  taken 
simuHaneously  with  the  O3  data  of  Figure  7.  The  small  quantity  of  emission  at  the  lowest 
SZA’s  is  from  a  blend  of  OH(AV  =  l)  and  CO2  2.7  /xm  bands;  the  increase  starting  around 
94°  is  ascribed  entirely  to  OH(V)  formed  from  the  increasing  quantity  of  ozone.  The 
SHARC  calculation  yields  slightly  more  OH  emission  than  observed,  but  the  difference  is 
well  within  the  variability  of  the  OH  emissions  and  the  pointing  uncertainty.  The  solar 
angle  dependence  of  the  data  is  also  reproduced  quite  well;  the  location  of  the  terminator 
edge  appears  to  be  accurate  to  within  2°  or  less  of  SZA. 

Since  the  concentration  of  ozone  derived  from  the  photochemical  model  is  inversely 
proportional  to  the  H  atom  concentration,  the  quantity  of  OH  emission  at  night  does  not 
depend  on  the  assumed  H  atom  profile.  However,  it  does  depend  on  the  O  atom 
recombination  rate,  and  hence  on  the  O  atom  concentration  as  well  as  the  mechanism  and 
rate  constants  for  OH(V)  deactivation.  Future  refinements  of  the  OK(V)  kinetics  in 
SHARC  will  provide  a  more  precise  relationship  between  the  OH  emissions  and  the 
recombination  rate,  so  that  these  emissions  could  be  used  for  deriving  accurate  O  atom 
profiles  in  the  mesopause  region. 
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SZA  (degrees) 


Calculated  (SHARC-3)  and  Observed  (CIRRIS-IA)  OH(AV=l)  Limb 
Radiances  at  the  Dusk  Terminator. 


4.  CONCLUSIONS 


The  SHARC-3  atmospheric  IR  radiance  code  has  been  used  to  simulate  limb 
radiances  for  CO2,  O3,  and  OH  in  the  solar  terminator  region  that  correspond  to 
observations  from  the  1977  SPIRE  and  1991  CIRRIS-IA  experiments.  These  calculations 
represent  the  first  direct  comparisons  with  field  data  using  the  new  multiple-profiles  and 
multiple-atmospheres  features  in  SHARC-3.  The  good  agreement  between  the  calculations 
and  the  data  provide  substantial  validation  of  the  kinetic  and  radiation  transport  models  in 
SHARC  as  well  as  the  O3  and  H  atom  profiles  from  the  SHARC  Atmosphere  Generator 
(SAG). 

Similar  comparisons  are  underway  using  additional  data  from  CIRRIS-IA  and/or 
other  exneiments  to  more  fully  investigate  these  emissions  over  a  wider  range  of  viewing 
conditions  and  to  investigate  other  radiators  and  bandpasses  of  interest. 
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